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I. INTRODUCTION 

The properties of multi quark-antiquark interacting systems as investigated by Weinstein and Isgur [l[ pointed to a 
possible bound kaon-antikaon molecular state structure for the light scalar mesons /o(980) and ao(980). Alternative 
proposals include, in addition to the molecular picture [U, H, H, 0, B @1 j a qq state Q or a q^q"^ state [1]. However, the 
large 27r annihilation width for qq ^ nn ^ 500 MeV from flux tube-breaking to ~ 600 MeV using current algebra 
[lOj, disqualifies a light quark-antiquark configuration when compared with a typical experimental width of ^ 50 
MeV for the /n (980) mr decay [lH. The inclusion of strange quarks changes this picture For example, it 

is shown in [I4I that if the /o(980) is identified with the (almost) pure ss member of the nonet within a quark-level 
linear a model framework, one can obtain reasonable values for the mass as well as tttt and 77 decay widths. So the 
question of a molecular versus a qq elementary particle structure, or perhaps some combination of these, for the light 
scalars remains an open one [3| • 

In the following we re-investigate the molecular state option using the non-local extension of the vector meson 
exchange potentials derived in [5| for kaon-antikaon bound states and their decay modes. It is shown that if one 
adheres strictly to the approximations required [l^ [11] for reducing the instantaneous version of the Bethe-Salpeter 
equation to a non-relativistic wave equation, especially as regards the off-shell nature of the relative four-momentum 
of the boson pair undergoing binding, a quantitatively similar picture to that given by a purely local approximation 
to the potential emerges for kaon-antikaon bound states. In particular, the non-binding states reported in a recent 
analysis [l3] of the same problem now no longer occur. 

The paper is arranged as follows: Section II discusses possible kaon-antikaon bound states of the coordinate 
space version of the non-relativistic Bethe-Salpeter equation for non-local interactions as generated by vector meson 
exchange from a S'[/y(3) x S'[/^(3) invariant interaction Lagrangian. In Section III we calculate the dominant decay 
widths of such bound KK pairs of good isospin for annihilation into either tttt or Tr^yy, while in Section IV the low 
energy KK scattering parameters in the presence of annihilation are found. Finally, Section V is devoted to a brief 
discussion of the energy shifts and decay widths to be expected for the associated mesonic atom K'^ K~ (kaonium) 
in context of the non-local interactions introduced in this paper. A summary and conclusions follow in Section VI. 



II. KK BOUND STATES 



A. Non-relativistic Bethe-Salpeter equation 

As before |5| we use the non-relativistic version of the Bethe-Salpeter (BS) equation in the instantaneous approx- 
imation [11,113, described further below, to study the mass and decay widths of the /o(980) and ao(980) scalar 
mesons, considered as weakly bound kaon-antikaon pairs of isospin zero or one interacting via vector meson exchange. 
In momentum space this equation reads 

where 0(p) is the wave function of the interacting pair in momentum space; Mk ~ 496 MeV is the average kaon 
mass. The eigenvalue Pq = 2Mk + E is the as yet unknown total mass of the interacting system of binding energy 
E^-e<0. 

Apart from introducing non-relativistic particle propagators, the other essential simplification used in deriving 
Eq. ([T]) is the assumption that the interaction is instantaneous. This is implemented by suppressing the time com- 
ponents of the relative four-momenta p — {pi — P2)/2 [0,p] and p' = {p[ — P2)/2 [OiP'l in the irreducible 
four-point vertex, or transition amplitude (pi,P2|r|pi,p2) of the original BS equation [isj . For a bound state the 
four-momenta that enter and leave this vertex are all off their mass shell. In the instantaneous approximation in 
particular this means that {pi,p2) — [|^o,±p] and (^1,^2) = [■^^Oiip'] in the center of mass (cm.) system, where 
Pi + P2 = Pi + P2 = [^OjO] still guarantees total energy-momentum conservation for arbitrary values of p and p'. 
These approximations result in the transition amplitude 

{p[,p'^\t\pi,P2) = (p',Po|fb,Po) «f(p',p;Po) (2) 

which appears in Eq. ([T]) that only depends on the unconstrained three-momentum of the incoming and outgoing 
kaons in addition to the total energy, and legitimizes integration over these independent momentum variables. 

An expression for f (p',p;Po) for KK scattering via t-channel vector meson exchange has been derived in [1, Q 
and also fl7j from a standard S'?7y(3) x 51/^(3) invariant interaction Lagrangian [18] to which we refer for further 
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details. One finds in this case that {p'i,p'2\i'\pi,P2) is given to lowest order coupling by 

-(Pl +K) • iP2+p'2) + {Pl-P?)ipl-P?)/Mf 



that reduces to the approximate transition amplitude introduced in Eq. 



(3) 



{p'i,p'2\^i\pi,P2) ~ fi(p',p;Po) 



+ (P ' + P)' + (P 



. 2\2 



(P ' - 



(4) 



for the values of the off-shell four-momenta prescribed above. This vertex structure also holds true for KK — > tttt 
where p' now refers to the outgoing pion momentum, but not for KK n'^r] without modification, see below Eq. (j26D . 

The coupling constants gf are all related to the the pmr coupling constant gpTnr by the SU{3) symmetry 18]: 
9t = {k,hh)9%7r for « = {P,^,(l>) exchange in KK KK, or gf = {glKK-^ 9^KK'gvKK') = (^,^V3)g2^^ for 
K* exchange leading to KK tttt or tt^t] transitions. Mi is the mass of the exchanged meson in each case, and 
as = gp^T^/in « 2.9 is known from the KSRF relation The Ci = [3, 1, 1] or [—1, 1, 1] are isospin factors [sj for 

(p, u), 4>) exchange between a {KK)i pair coupled to / = or 1, or d = [— \/6, — \/2] for the / = and 1 transitions 
[KK)f)_i — > (7r7r)o or (7177)1 via K* exchange; {KK)i — > (tttt)! is forbidden by G parity conservation; Us is a symmetry 
factor [-^, 1] for identical (non-identical) bosons in the final state. 

As already remarked in , one can check the relations between the various vector-pseudoscalar coupling constants 
given above empirically by appealing to the experimental decay widths of those vector mesons i = {p, K* , cp) that are 
unstable with respect to decay into tttt, Ktt and KK. An elementary calculation shows that a vector meson of mass 
Mi at rest has a decay width of 



3M? Att 



(5) 



in the scheme of Q to first order in the relevant coupling constant gf, where pi is the common cm. momentum 
magnitude of either decay particle in the final state, and fi = (2,3,2) an isos pin factor for the three channels in 
question. The symmetry factor Ug has been defined above. Inserting the known masses and full or partial decay 



widths of {p, K* , (p) as appropriate into Eq. ^ , one finds {g 
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9. 



4>KK' 



(1 : (0.54)2 . (0.77)2) compared 



with (1 : (1/2)2 . (^.^2)2) given above, where gp^^/An = 2.92 that coincides almost exactly with the KSRF value. 
The remaining empirical ratios hold to within ~ 8% of the "ideal" mixing version of the SU{3) relations between the 
coupling constants. 

The expression for the two-particle vertex given by Eq. ([4]) regards the basic vertices as point-like. Upon introducing 
a form-factor depending on the three-momentum transfer q = p' — p of the form 



, /2A2-A/2n2 



(6) 



at each of these vertices with a common high-momentum cutoff ^ V2A, the vertex in Eq. ([T]) is altered to read 
r(p', p; Pq) = J2i ri(p', p; Po)Fi{l^)- The use of a common A is simply a working assumption; the cutoffs could in 
principle be different for different vertices [20| . 



B. Non-local interaction potential in coordinate space 

The vertex in Eq. ^ was approximated further in ^ by neglecting the small binding energy e ~ ^-^if contribution 
by taking Pq « 2Mk and ignoring the three-momentum dependence in the numerator entirely. Then the components 
of the KK potential on the right hand side of Eq. ([T]) reduce to 

-l^t,ip',p;Po)Ff{ci^)^gfQU.iq), U,i<i) ^ (7) 

that is only a function of the momentum transfer q'^ . This leads to the following set of local potentials in coordinate 
space for each exchanged meson Mi that are known in closed form from the fourier transform of Ui{q), 

1 —Mr — \/2Ar 

V.ir)^-agfU.{r), [/,(r) = -{^ - J] C(°)(y2Ar)"}. (8) 

n=0 
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The d°^'s are polynomials C^°) = 1, Cf^ = j^(ll-4cf + c4)(l-c2), = ^(3-c2)(l-cf )2 and ^ = i^{l-^,f 
in the variable Cj = Mi/y/2K. The partial potential Vi{r) has a simple power series expansion in r about the origin 
when the form factor is included, 



V,{r) = -V2Aa,ail - c.)'{^(5 + 4c, + c^) - 1(1 + 4c, + c^){V2Arf + • • • } 



(9) 



with ai = gf/4:n, so that the Vi{r) as well as its space derivatives are all well-behaved at r = for finite A. 

In general the interaction vertex in Eq. ([T]) leads to a non-local potential in coordinate space. However, since the 
numerator in the special case of Eq. ([5]) is a polynomial in the three-momentum variables, one can incorporate them 
in coordinate space in a revised potential that also involves space derivatives of Vi{r) coming from the replacements 
p or q ^ — iV in ri(p + q, p; Pq)- Again setting Pq « 2Mk, the KK non-local interaction potential operator V^r 
becomes 



1 



AMI 



V2{r)V^ + 2V • V2{r)V + VV2(r) 



(10) 



when operating on the s-state wave function tA(^)i the three-dimensional fourier transform of the spherically symmetric 
momentum space wave function 4i{p) of Eq. ([1]) . 

The potentials V'i{r),V2{r),Vi{r) are given by the following combinations (the primes indicate derivatives with 
respect to r) 

1/3 (r) = Fi(r) + ^i(F/(r) - Vi{r)) 



V2{r) = 

i 

Vi{r) =Y.V,{t). 



Vl'jr) 
M? 



(11) 



Then the coordinate space version of Eq. ([T]) with the full momentum structure of Eq. ([4]) included can be written as 



1 



■2V- 



1 



v + 



.M^(r) M^(r) 
containing a spatially-dependent "effective" kaon mass 



M^(r) 



M 



^-Kir)=l~\r)MK, 7(r)-(l + -^) 



(12) 



(13) 



that enters into the equation via a fully symmetrized kinetic energy operator. 

The quartic terms in the vertex Eq. ^ only contribute to off-shell scattering; they vanish on-shell when p'^ = p^. 
This suggests that their contribution to the interaction potential may be small. If they are omitted entirely, both 
V3 (r) and V2 (r) reduce to Vi (r) and the structure of Eq. becomes identical in form with that of a wave equation 
derived in ^21] that has been used quite generally for describing non-local effects on particle motion in nuclei [2^ . 

The hermitian structure of the symmetrized kinetic energy operator assures that the eigenstates of Eq. remain 
orthogonal, and also that the continuity equation for its time-dependent version continues to hold with a revised prob- 
ability current density j(r, t) = \ / iM'^{r)[ip*'S/ — ij)*]. Thus the normalization ofip(r,t) also remains independent 
of time. 



C. Numerical results 



We next investigate possible bound s-state solutions of Eq. p^ . It is convenient to set ip{r) ~ ip{0)u{r)/r with 
u{r)/r ^ 1 as r ^ 0. The value of ■(/'(O) is fixed by normalization. Then Eq. p2)) simplifies to 

7(r)w" + i{r)u' + Mk{E - Vo)u = 0, (14) 

where 

y^ir)^V.ir)^^\viir)-^Vi'ir). (15) 
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We solve this equation numerically after using the partial potentials given by Eq. ^ to construct the required 
interactions. For / = one has 



and 



Vi{r) = ^g^.QUpir) + lc/^(r) + ^U^r) 



M,2 



Mi 



(16) 
(17) 



from which Vb(r) can be obtained. For 1=1 identical expressions hold with the coefficient 3/4 of the p meson 
exchange contribution replaced by —1/4. 

Since known physical masses of exchanged mesons {Mp, M^, M^) = (768,783, 1019) MeV and their coupling con- 
stants given above enter the calculations, the only free parameter is the cutoff A. The possible choices of A are 
then restricted further if one in addition requires that 7(0) > 0, in order to ensure that M^(r) remains positive in 
the interaction zone [s^l- This problem has already been encountered in [l3l, where it is concluded that for cutoffs 
that ensure (their equivalent of) 7(0) > give non-local potentials that do not bind the KK pair. However these 
potentials have been based on an approximated vertex that corresponds to ri(p', p; Po) of Eq. ^ without the quartic 
contribution and Pg replaced by the combination 2{Ep, + Ep) of kaon cm. energies. The first approximation neglects 
a small correction term, but the second one excludes off-shell kaons, p' p, hy total energy conservation, contrary 
to what is expected for bound states. Thus these potentials may not correctly reflect the underlying physics for such 
systems (isj . 

This is not the case for the potential operator Vj^j^ of Eq. (|10p . Using Eqs. Q and (llip one readily finds that 



7(0) 



1 



V2A 



J2 c^^C.i^ - c0'{^(5 + 4q + c\) + ^(1 + 4c7i + cr^)}. 



(18) 



For the known interaction strength of = g^^^/47r = 2.9, the factor 7(0) stays positive in both isospin channels if 
A lies in the interval 390 ^ A <^ 1280 MeV, see Fig. [TJ However whether or not the resulting non-local potentials 
bind the KK pair is still sensitive to the actual choice of cutoff. Table Uhsts calculated values of binding energies, 
masses and decay half-widths (see next section for the latter) for the subset of cutoff values A = (415, 420, 425) MeV 
confined to a relatively small window within this interval. These cutoffs have been chosen so as to bracket the spread 
in experimental [ll[ binding energies 11.3 ± 10 MeV for the /o(980) meson considered as an interacting KK pair. As 
a byproduct the ao(980) then also appears as a bound pair state in the isovector channel of these potentials. 

A much stronger purely local potential Q of depth ^ 4 GeV was required to give the same order of binding 
for / = (and none at all for / = 1). The reason for this is clear from Table [TTl which summarizes results for 
the isoscalar potential. We see that the much shallower non-local potential ^ 0.2 GeV is counter-balanced by a 
significant suppression of the kinetic energy through an increased effective kaon mass of M^(0) « 2Mk at the origin 



that favors binding. One notes that cutoff momenta 



\/2A ~ M,,, of order of the 77 meson mass, that also sets 



the mass scale of the pseudo-scalar meson octet, are required to regulate the behavior of the bare vertex in Eq. ^ 
at large momentum transfer. Otherwise this vertex either tends to a constant or diverges depending on whether 
the fourth order contribution is dropped or kept. These cutoffs are typically an order of magnitude smaller than 
those for a local potential like that in Eq. ([S]) to give the same order of binding. We also comment that the results 
in Table H] remain qualitatively unchanged by dropping the quartic momentum contributions altogether and setting 
V^{r) = V2{r) = yi(r) in Eqs. (HH) and (US]). 



III. DECAY OF KK BOUND STATES 

The annihilation of bound KK pairs into mesons or photons has been discussed in (5| and [2^, [131 in the context 
of the BS equation. The decay of the bound state is governed by the imaginary part of the additional box diagram 
rf,ox(p', P; Po) contribution to the BS equation vertex where a K* vector meson is exchanged in the t-channel, see 
Fig. [2I Then the eigenvalue Pg in Eq. ([T]) acquires an imaginary part, Pn — > Po — iV /2. The imaginary part of the box 
diagram can be retrieved directly by applying the Cutkosky rules [l^ to cut the intermediate meson loops in Fig. O 
The decay widths F are then obtained perturbatively by averaging over the bound state solutions of Eq. ([1]) in the 
absence of annihilation. 

For the / = 0, KK tttt decay this width reads 

^- = 32^^44/'"^'"'"^"'^^'' ^''^ 
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where the bound to free transition ampHtude MsipTr) can be written as either a momentum or a coordinate space 
integral, 



(20) 



by fourier transforming 0(p) and MttttIPtt, p)- The angular integral drijr in Eq. runs over the full solid angle of 
one of the pions due to wave function symmetrization [3j] . 

The amplitude AfjrTrlPTr, p) refers to KK annihilation where the pions are on-shell with the magnitude of their 
common momentum fixed by total energy conservation at p^r = {jPq — w^)^^^; the kaons remain off-shell as before. 

M^7r(P7r,p) is then given by f i(p^,p;Po) + f ,j(-p^, p; Pq) from Eq. (g]) with gfCiUs = -V^glxK- and Mi = Mr-- 
The p^ — > — Ptt crossed contribution comes about since the pions are identical bosons in the symmetrized space x 
isospin basis |3[ . Performing the fourier transform of this sum gives the transition operator in coordinate space as 

M^AP.,r) = -2V3glKK'e~''''-'\sHr) + [P^ + ^pI - Mj,, + 4z(p, • V)J ^— j (21) 

after discarding the quartic terms of order C'(Af^/4M^,) ^ 1. The delta function contribution arises from the 
quadratic dependence on the momentum transfer q = p — p^ in the numerator of Eq. ^ that introduces the 
Laplacian operating on the Yukawa-like potential in coordinate space produced by the K* meson exchange. We 
insert this result into the second form of Eq. ([20)1 with ijj^r) = ^{0)u{r)/r to find 



MeiPn) = -87ra^^V(0)i?(p^) 
where a^^r is an effective coupling constant at zero kaon momentum 

r Po+pI 



r^(p.,0;Po) 1 /TT 

= -V3 as 

47r 4 



Mf,.+Pl 



1.489 



(22) 



(23) 



for Us — 2.9 and Pq ~ 2Mk- This determines the annihilation cross section for free KK pairs into two pions at low 
momentum as 



Ml p- 



The factor R is given by 



R{p.) 



Mh+pl 



l + iP,'+4pi-Mf,, 



dru{r)jo{p^r)e 



^pI 



The ji{p-Kr) are spherical Bessel functions. Then Eq. (|19p reads 



(24) 



'dru(r)(l + MK.r)^i^^e-*^^"^')j (25) 
p^r )\ 



Ml Mk 



)V'2(0)p2(p^) = {vM)i^''{0)R'{p-.) 



(26) 



where Vrei — "^p/Mk is the relative cm. velocity of either kaon and Pt^ = 0.959 Mk is now the pion momentum at 
threshold. For P = 1, Eq. ((26)) reduces to the familiar "wave function at contact" approximation for the decay width 
[2,, 5, 26, 28i] where only the probability per unit volume '0'^(O) of finding the pair at the origin characterizes the bound 
state. 

The isovector KK — > Tr'^ry decay width F^o^ is given by a similar expression with the following replacements in 
Eqs. (121 to (ESD: p^ by p^^ = 0.658 Mk, P^ by (P^ + iPo)(P^o + iPg), and Mj., by M^, = M]^, + (P^ - 
iPo)(p7rO — \Po)- These modifications of the vertex in Eq. (H reflect the revised total energy conservation condition, 
Pq — P^o + P^ for an on-shell -K^rj pair in the final state with energies E^^o and and common cm. momentum 
Pi^ri- Then a^rTr is replaced by 



{Er, + iPo)(P.O + iPo) 



Ml, 



1.079 



(27) 
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after removing the crossed contribution; 'ip{0) now refers to the isovector channel. 

In Eq. ipS]) both ■0(0) and R are functions of the binding energy; R in particular only differs significantly from unity 
when the ranges of the radial wave function u{r) and the amplitude M^rTr (Ptt , r) under the integrals in Eq. ([25|) are 
commensurate. This is the case for the analogous Mjj{p^,r) that describes two-photon decay [26|. However for the 
loosely bound KK pairs in Tables HI and HH u{r) has a typical range ^ that is about twice the range M^l of the 
transition operator. Then only the short range behavior of u(r) is important for determining R. This behavior does 
not differ very much from the small r limit, u{r) = r, when both integrals have the common value (Af^, +p^)~^. This 
gives R — 1, thereby reproducing the wave function at contact approximation. The same result follows immediately 
from the momentum space representation of Mb(j>t:) in the weak binding limit when (/)(p) has a much shorter range 
than the transition operator in momentum space. Then the integral evaluates approximately as Mt;^{p„, 0)0(0) that 
leads directly back to Eq. ([^^ again with R—1. 

The calculated widths are listed in Table [J Since the cross section factors Vrei^'i for free annihilation are fixed 
by the coupling constants a-^T^, ctTvrj, the values of the decay widths given by Eq. ^2E\i and its / = 1 counterpart are 
controlled by the real KK potential that determines '0(0) and R in each isospin channel. The cutoffs A have been 
fixed to reproduce a prescribed binding energy range. There are thus no free parameters. The predicted widths in 
Table I cover a range for F^^ that agrees quite well with experiment within the quoted error bars; the predicted values 
for r^o,^ are too small. 

The actual values of R listed in Table HIl for / = increase only slightly with binding energy from 1.034 to 1.041 over 
the entire range of cq. The behavior for / = 1 is similar. Thus the main variation in width in Table U comes from the 
contact probability density 0^(0) which increases approximately like higher binding enhances the probability of 
contact at the origin. From a quark model perspective this increased contact probability facilitates strange quark 
annihilation and exchange to produce the {KK)q tttt and {KK)i tt^t] decays respectively. 



Contributions from s-channel scalar meson exchange 

The calculations of mass and widths just presented for the KK system ignore any s-channel contributions to the 
scattering amplitude r(p', p; Pq) in Eq. ([1]). However, KK can also interact and annihilate via scalar meson exchange 
in this channel. One expects such effects to be small because the exchanged particle is off-shell. As an example we 
estimate their order of magnitude contribution for the case of a meson exchange, treated as a simpl e qq st ate. Then 
the relevant interactions can be read off from the SU{3) linear sigma model (LcrM) Lagrangian [12. [isL l29l|. 

Ccubic = gaTnrO-{n ■ tt) + g^j^f^aKK H (28) 

The couphng constants are given by ga-^T: ~ M^/2f.^,g^j^j^ w M^/2fK in the (ideally mixed) chiral limit. The 
(/tt, Jk) ~ (92.2, 110) MeV with fK/f-n ~ 1-19 are the pion and kaon weak decay constants [TT'|. 

The relevant diagrams are shown in Fig. [3) Translating these diagrams one finds a s-channel contribution to the 
BS equation scattering amplitude coming from a exchange as 

f4p',p;Po) = -{V2g,KKf'DAPS)) (29) 

for off-shell kaons with cm. four-momentum squared s = Pq , where I?ct(s) = — — Ecr(s)]~^ is the "dressed" 
scalar meson propagator [l^ of mass and self-energy Eo-(s). 

The s-channel contribution to the additional binding energy and decay width for / = is then given in first order 

by 

In particular, 

4';^=^%^^'(OW(^o), pA^) = ~-ImV,{s) (31) 

where Pa{s) is the scalar spectral density associated with 'D„{s). 

The self-energy Yjc{s) receives contributions at one-loop level from tadpole and seagull diagrams in addition to 
polarisation diagrams generated by Eq. (j28p . plus counter-terms. Of these, only the polarization loop involving vrvr for 
/ = (or TTjy for / = 1) has an imaginary part for strong decay at the relevant four-momentum transfer encountered 
in Eqs. (|30p and (|3ip . Cutting the loop in the standard way places both pions on-shell and gives the imaginary part 

- /ml],(s) = 2^(1 e{s - AMI) = M,r,(s), (32) 

ovr V s / 
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The last equality defines an "off-shell" width ro-(s) that coincides with the physical decay width for cr tttt at 
s = M^. Evaluating p^is) in Eq. ([?T|) at s = Pq, one has 

^f^) ^ 3p^^'(0) f 91..91kK \ /onN 

after absorbing i?eScr(s) into a redefinition of the a mass AfJ + (i?eEcr(^'o ))^ ^ A'^^ • The cm. momentum refers 
to outgoing pious from the KK tttt decay as before. The binding energy contribution e'^^ follows from the real 



part of Eq. ^ as 



{P^ Mj)[2/m5].(Po2)]-ir(';) < 0, (34) 



i.e. the contribution to the real potential from the s-channel is repulsive (less binding) for Pq > M„. 

We identify a with the /o(600) scalar meson and take its mass and decay half-width parameters as M„ — — 
[(541±39)- (252±42)i] MeV directly from experiment By inverting Eq. ^ for the a ^ nn width the couphng 
constant is determined as go-n-Tr = 1-63 ± 0.17 GeV. Since this value is quite close to the chiral model estimate of 
ga-KTT = Ma/2fTr ~ 1.48 GcV, wc assumc that the LtrM coupling constant ratio {gaTn^)/{gaKK) = Ik/ fir ~ 1-19, also 
continues to hold to fix g^^R — 1-37 GeV. These parameters lead to / = contributions from the s-channel of 

4'"^ + ^r^;;' = -8.92 + 3.93i, -5.18 + 2.28i and - 2.92 + 1.29i MeV (35) 

respectively for cutoffs A = 415, 420 and 425 MeV. 

The exchange of other scalar mesons also has to be addressed. Since both /o(980) and ao(980) are depicted as KK 
molecules in the Weinstein-Isgur picture used here, the next available set of scalars to consider would be the O"*" states 
/o(1370), iirg(1430), ao(1450) and /o(1500) with masses above 1.3 GeV. However a consistent theoretical description 
of these scalars is still under debate. If, as suggested in [H, [H,!!^, one places them in the same qq scalar nonet and 
implements a linear sigma model description, the decay widths of, for example, ao(1450) into KK and tt'^i] come out 
far too large. This is in part due to the large coupling constants generated by the model [13]. On the other hand if one 
assumes an interaction of the LctM form, but extracts the relevant coupling constants (gaoTrri, gaoKii) ~ (1.34,0.95) 
GeV from the ao(1450) branching ratios iTlJl into tt^t] (^ 8%) and KK (^ 7%), the I — 1 contributions are given by 

e!""'' + ^r^i"J = 2.16 + 0.085i, 1.15 + 0.045? and 0.59 + 0.023i MeV (36) 

for s-channel ao(1450) exchange. We omit the calculational details. 

None of these contributions introduce any significant corrections into the t-channel values of Table HI This is because 
the scalar spectral density is sampled at s « 4M|- that pushes the exchanged meson significantly off its mass-shell 
to weaken the s-channel KK interaction accordingly. In addition, the small partial decay width ^ 20 MeV for 
ao(1450) — > ir'^r] relative to u — > tttt serves to suppress the 1=1 contribution even further. 

In closing this section we remark that the interaction vertices described by the effective Lagrangian involving only 
meson degrees of freedom can also be visualised at the constituent quark level. For instance the vertex if + + K'^ * — > 
7r+, that contributes to K~^K~ — tt+tt" in Fig. [2] translates schematically into [{us){sd)] the strange quarks 

ss having annhilated (via some unspecified gluon exchange interaction) leaving behind the u and d quarks. 

Moving beyond the effective meson theory model then, this suggests that in principle the final two-meson decay 
channels can for example also be reached via another route involving "constituent gluons" in intermediate states that 
give rise to hybrid meson structures [qqg] (see [3]| for a comprehensive review). Hybrids are generally expected [3^[33| 
to be heavier than 1.7— 1.9 GeV, thus lifting the threshold of such intermediate states above 2 GeV and making them 
less important than the intermediate meson states which we have taken into account in the present context. Moreover, 
the hybrids' predicted decay properties (for example no decays into two pseudoscalar mesons, at least according to 
the gluon flux-tube model [34S of their structure) can only further reduce their impact here. Further discussion of 
such contributions lies outside the scope of this article. We refer to the recent literature [U, [12, [13, [IJ for more 
details. 



IV. KK SCATTERING 

It is straightforward to show that the standard effective range expansion [1^ for the KK scattering phase shift 
Sj{k) at cm. momentum k = ^/M^E, 

kcotSiik) = -— + -ri + ■ ■■ (37) 
ai 2 
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continues to hold for the modified radial equation ([14]) in terms of the scattering length and effective range aj and rj. 

In the presence of annihilation both parameters pick up imaginary contributions. In particular, in the limit fc — > 0, 
the imaginary part of the 7 = scattering length qq is related to the KK — > tttt annihilation cross section (Tq by 

-Im(a„) = |ja„ « A|/(0)|-V„ = 1/(0)1-^-' 

The first step is exact. The next step follows after notin g th at ctq is related to the free annihilation cross section ctq 
of Eq. p4)) by (Tq « |/(0)|"^(To to a good approximation [3^. Here f{k) is the Jost function [s^IstI t hat determines 
the scattering matrix S{k) = f{k)/f{—k). The value |/(0)|~^ at fc = is the enhancement factor that gives the 
ratio of the probability of finding an interacting kaon pair at the origin to that when there is no interaction between 
them. 

If one parametrizes the Jost function /(fc) = |/(fc)| exp[iSi{k)] as 

fik) = ^ (39) 
k ~ lb 

the effective range expansion becomes exact [37j . and scattering length and effective range are given in terms of the 
parameters (a, 6). Considering 7 = again, one gets 

ao^-{b~a)/ba, ro = 2/{b-a). (40) 

Since the isoscalar channel supports a single s-wave bound state at complex binding energy sq — sq + ^Ftttt one 



knows that the function f{k) must have a zero 37| in the lower half of the complex k plane at fc ia — —i^/MKCo 
which determines a. Combining Eqs. ([55]) and (PD|) one then finds b — \J ij^Im{—a) for b real. An identical procedure 
suffices for determining (a, 6) in the isovector channel. 

The scattering parameters for both channels are then given by Eq. (|40p . Taking the relevant values of eg and ei 
from Table I for A = 420 MeV, one calculates that 

ao = 4.926- 2.235J, tq = 2.550 - 0.525i 

ai = 8.246- 2.576J, n = 2.971 - 0.242i (41) 

-1 



in units of Af^ 



Apart from one pioneering experimental attempt to determine the isoscalar scattering length by Wetzel et al. 
who find ao = [|(3.13 ± 0.30) | — (0.67 ± 0.07)i]Af^^, and a subsequent analysis of later tttt data that infers 
ao = (4.36 — 1.49i) M^^, there are as yet no other direct KK measurements with which to compare the estimates 
in Eq. (|41[) . However, to the extent that only coupling between the tttt and KK channels is important, two channel 
unitarity shows (40| that the tttt and KK isoscalar inelasticities rjo^mr) — rio{KK) — |5o(fc)| are equal at their common 
total cm. energy Pq = 2-\/p2 + Af^ = 2-\/fc2 + M"^ where p and k are the pion and kaon momenta. In Fig. 3] we 
compare rjo^i:?:) as calculated from S'o(fc) — exp[2j(5o(fc)] in the effective range approximation, Eq. (j37[) . using the 
isoscalar parameters given above, with the inelasticities extracted from tttt scattering in Ref. [41]. The result using 
the isoscalar scattering parameters ao ~ (4.281 — 2.398i)M^^ , vq — (1.169 — 0.178i)Af^^ obtained from the local 
potential in Q is also given for comparison. A recent TC-matrix fit based on a combination of various data sets is 
also shown, see (42j and further references cited therein. Given the very wide error bars on the data it is difficult to 
draw definitive conclusions beyond remarking that the model calculations do reproduce the correct trend and order 
of magnitude of the tttt inelasticity, albeit somewhat too low. 



V. KAONIUM 

We next summarize the results obtained for the energy shifts and decay widths for the K~^K~ (kaonium) atom in 
the context of the present calculations. The properties of this system have already been discussed in detail in Q for 
local potentials. 

Pure Coulomb interactions bind kaonium at — ^a^/u = —6.576 keV in the lowest Is state where V'is(^) = 
7r~^/^(^Q;)'^/^exp(— ^ar); fj, = ^Mk± is the reduced mass with Mx± ~ 494 MeV, and a ~ 1/137 the fine struc- 
ture constant. The kaonium — > tttt + n^rj decay width follows from Eq. ([26]) after replacing ctq by the total K^K^ 
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annihilation cross section ap = (cto + ci)/2 and ■0(0) by i/jis^O). While R can be calculated in closed form for the Is 
Coulomb ground state of kaonium, 





-pI 1 


-1 r 









+pI- fia{2MK' + 3^q) ^ ^fia / {Mk- + mY 

{Mk' + fia^ +pI Ptt ^ Pt, ' 



(42) 



this factor differs but little from unity, Ris{pt:) — 1.0016; the corresponding expression for 1 = 1 gives RisiPnri) = 
1.0013. This is to be expected. The Bohr radius of kaonium, 1/iia = 109 fm, is about 500 times the transition 
amplitude range of M^l - 0.2 fm. Hence the wave function at contact is a very reliable approximation in this case 
too and yields a total Is decay width of [FH 



Tis = + r^o,, ^{al^p^ + al^p^„)a^ = (51.3 + 18.4) cV = 69.7 eV, (43) 
or a lifetime of ^ 10^ ^^s. 

The effect of strong plus Coulomb interactions on the kaonium spectrum can be included in the standard way (35l | 
by equating the logarithmic derivative of the asymptotic form Up(r) ^ 1 — a^r of the zero energy scattering 
wave function outside the strong interaction zone with that of a pure incoming Coulomb wave function; ap is given 
below. Physically this suggests that the binding of kaonium is still essentially provided by the long range Coulomb 
attraction. 

For bound states the Coulomb wave function is proportional to the Whittaker function [43| fi~\k,r) ~ 
Wij^ i/2{'2ikr) at complex wave number k\ = —iX^a that describes a decaying wave at r —^ oo. Here A is a yet 
to be determined eigenvalue that gives the binding energies and total decay widths E\ — = —^X^fia^ of the 

mesonic atom under the combined influence of strong and Coulomb interactions; irj = fia/ik\ = 1/A is the Coulomb 
parameter for attractive interactions. The matching condition at the Bethe "Coulomb joining radius" (ssj d say, is 
then equivalent to the Kudryavtsev-Popov equation [i^l for the complex eigenvalues A that takes the simple form Q 

{ap)c - 2fia [V^l - 1/A) + ^A + In A + 7] (44) 

for kaonium in the limit ^ad « 1, where ip is the digamma function [i^ and 7 = 0.5771 • • • is Euler's constant. The 
{ap)c on the left is the Coulomb corrected K~^K~ inverse scattering length |3a |. 

(ap)c — ctp — 2^a[hi{2^ad) + 7]. (45) 

Here ap = l/op is the inverse of the K~^K~ strong scattering length ap = (ag + ai)/2 without Coulomb corrections 
if the isospin-breaking arising from the kaon mass difference A = Af/^o — M^± w 4 McV [lH is ignored. If not, the 
K^K^ K^K^ charge exchange channel is closed, and [i^ 

/^-fcoa^y fco = v/2M^. (46) 
\ flp — KQaoai / 

We take the scattering lengths of good isospin from Eq. ([1T|) . and calculate the Bethe joining radius as d = 1.636M^^ 
from the regular zero energy scattering solution of Eq. p4)) . Then ap = (6.873 — 3.387i)~^MK and (ap)c = (5.948 — 
2.364i)~^Mif . Inserting this information on the left hand side of Eq. (j44|) one finds Xis = 0.9579 + 0.0161«, or 

Els - ^^is = (-6.027 - 0.205i) keV (47) 

for the lowest state. Repeating the exercise for the remaining two values of the cutoff, the results for the energy shifts 
and decay widths for kaonium ground state can be summarized as follows. 



AEis - ^Tis ^ {QMQto^ill - 0.205i[3:J?|ji) keV, A = 420 ± 5 MeV (48) 



that illustrates the sensitivity to the choice of A. The corresponding lifetimes read r = 1.6^°'qq x 10~^^s. 

Thus for A = 420 MeV, apart from causing a repulsive level shift of AEis = 0.549 keV in the ground state energy 
of Coulombic kaonium, the strong interaction also enhances the decay width considerably to Fi^ = 0.410 keV over 
that calculated in Eq. There is thus strong mixing between the pure Coulomb state and the /o(980) molecular 

ground state that produces a significant increase in decay width (and energy shift of the same order) that shortens 
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the lifetime of the kaonium ground state to 2 x 10~'^^s. These results are in line with previous estimates based on 
a local potential description [5|] . 

The strong mixing feature Fis ~ AE'is persists for the excited states of kaonium as well. This is illustrated for 
A = 420 MeV in Fig. [S] which shows the decay width versus the energy shift for the kaonium spectrum. The ratio 
is seen to remain remarkably constant over three orders of change in magnitude in these quantities as one moves 
through the Balmer spectrum of kaonium. This can be understood from the Deser et al. (47| perturbative solution 
of Eq. (I44p that is obtained by expanding the roots A^^ in powers of the corrected scattering length combination 
^a{a,p)c — ^a{ap)~^ around their pure Coulomb field values at = n = 1, 2, 3, • • •. Then 



— (ap)c[l + 2^a(ap)c(^VW - V'(l) + -Innj +•••] (49) 



that has a very similar structure to the analogous chiral perturbation theory result for kaonic hydrogen ^48i] . Working 
to lowest order, one retrieves the Deser formula 

A£;„, -^r„,«?^(ap)e = (0.575 - 0.228i)keV, n = 1, (50) 

which gives n — 1 shift and width values that are 5 and 11% too large in comparison with Eq. (|48p . The relation 
Tns = —2Im{ap)c/Re{ap)cAEns — 0.79A£'„s that follows from the latter result gives the straight line in Fig.[5l 

The underlying reason for the mixing can be qualitatively understood by temporarily ignoring annihilation so that 
Qfp is real and positive. Then one sees from the asymptotic form of Up{r) that the kaonium ground and excited state 
wave functions of Fig. [5] all develop a common additional node at r = 1 /ap that renders them orthogonal to the 
ground state wave function of the bound KK pair. They are consequently all excited states in the combined strong 
plus Coulomb potentials, that suffer level repulsion with the kaon-antikaon ground state. 



VI. SUMMARY AND CONCLUSIONS 



It is shown that the KK momentum-dependent interaction due to vector meson exchange can be incorporated 
without further approximation into a non-relativistic Bethe-Salpeter wave equation in coordinate space containing, 
in addition to a local potential, a spatially dependent effective kaon mass embedded in a fully symmetized kinetic 
energy operator. 

A numerical analysis of the possible bound and scattering states of this modified equation is presented as a model 
for the properties of a molecular kaon-antikaon description of the light scalar mesons /o(980) and ao(980). We show 
that for a restricted range of cutoffs A, the only adjustable parameter in the calculations, one obtains bound state 
masses for both scalar mesons in reasonable accord with experiment. The accompanying tttt and tt^t] decay widths 
exhibit a substantial variability with A since they are sensitive to the contact probability density as determined by 
the ground state eigenfunction. The results can be displayed in round numbers as follows 

Mo - 1f^^ = [(988 to 971) - ^(36 to 108)] MeV; Mi - ^F^o^ = [(989 to 976) - ^(7 to 27)] MeV (51) 

as the cutoff A increases from 415 to 425 MeV. The predicted range for Ftttt overlaps nicely with the experimental 
spread of values (40 — 100) MeV; the decay widths T^Oj^ for molecular ao(980) are too small. We also demonstrated 
at the end of Section III that the associated s-channel contributions to the binding energy and decay widths in both 
isospin channels are unimportant. 

These calculations treat the decay vertices as point-like. Introducing a form factor like that of Eq. © replaces i? ~ 1 
by Ra ~ 0.34 in the transition amplitude of Eq. ((22l) if we assume the same cutoffs again as used for the potentials, 
see Appendix Eq. (|A.7p . Then one obtains an order of magnitude reduction in widths down to F^j^ — (4 — 12) and 
FjrO^ ~ (1 ~ 3) MeV from those given in Eq. (I5T|) . This agrees qualitatively with the calculations in 49] that also find 
small decay widths of this order from a numerical study of the relativistic Bethe-Salpeter integral equation for the 
KK system when form factors are included. They in addition find that bound states only occur for cutoffs A lying in 
a series of narrow windows similar to the one illustrated in Fig. [T]for the non-relativistic case. 

The low energy parameters for KK scattering in the presence of annihilation have also been calculated in both 
isospin channels. In particular this information allows one to construct the KK isoscalar inelasticity as a function of 
the total energy in order to compare with the measured tttt inelasticity which it equals in the two-channel problem 
considered here. The qualitative agreement with experiment is satisfactory. 

Finally, the effect of strong interactions on the binding energy and decay of the K~^ K~ mesonic atom are briefly 
discussed in the context of momentum-dependent potential model. The energy shifts and decay widths so obtained 
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are again not very different from those using the local potential approximation, and give a lifetime of '-^ 2 x 10~"'^^s 
for kaonium. 

Overall, one can conclude that both the local and non-local potential description for the interaction of a KK 
pair lead to very similar results that are only dependent on the cutoff required in the form factor of the potential for 
obtaining bound states. However the width predictions of either of these, or other [i^] . calculations remain ambiguous, 
depending as they do on whether the vertices in the transition amplitudes are regarded as point-like or extended. 
The calculation from first principles of such hadronic size effects lies beyond the scope of the present work. Without 
a better understanding of this latter aspect, one cannot decide on the basis of width predictions alone between 
the molecular versus elementary particle interpretations of the dynamics governing the structure of the light scalar 
mesons. What has been shown here, however, is that the molecular model for /o(980) and ao(980) with either local 
or non-local meson exchange interaction potentials does give reasonable mass and width predictions for the isoscalar 
channel in particular, assuming point-like vertices in the decay amplitudes. In this regard, observation of the energy 
shift and decay width in kaonium, which in combination directly probe the strong K~^K~ scattering length, would 
be a particularly important source of experimental information for further progress towards a better understanding 
of these interesting systems. 
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APPENDIX 



Decay widths with form factors 

In order to include a form factor in the transition operator Mtttt (Ptt , r) of Eq. (PT]) we replace the pure Yukawa 
potential in that equation by Ui{r) from Eq. ^ with Mi = Mk' and q — Mk* / V^A. Then the revised transition 
amplitude splits into two parts. 



MA,..(p.,r)= -2V35^^^,.e-*P- 



Slir) - {P2 + 4pl + 4zp, • V}[/A(r) + {P^ + Apl - M|.. + 4zp, • V}- 



(A.l) 



after using 



V2[/,(r) = -Ml, 



-Mr,,r 



Airr 



(A.2) 



where 



(7a (r) 



-\/2Ar 

47rr 



^Cf(V2Ar)« 



n=0 



(A.3) 



and 



4(r) = (x/2A)2i^5:c(^)(y2Arr. 



(A.4) 



The coefficients are linear combinations of the ci°^ cf^ = {d^'^ -2Cf'^ +2d°^), cf ^ = (cf ^ - 4C^°^ + eC;^"^), 



while 



(6*2°'' — QC'^^) and Cg'^'' = C'^\ Since 6\{r) ^ or oo depending on the order of the limits A —f oo, 



^(2) _ ^(0) 



f dV5i(r) = ^C(2)r(n + 2) = l 

n=0 



(A.5) 



13 



independent of A, the function defined by Eq. (jA.4[) is a representation of a spread-out delta function for finite A. 
Clearly Ma^tttt (Pn- , r) reverts back to M7r7r(p7r,r) when A oo. For A oo Eq. (|A.ip shows that the factor R{pt^) in 
the bound to free transition amplitude of Eq. ([22|) is replaced by 





-pI 1 


-1 


[( 




+ pV 







3 3 »oo 

n=0 n=0 "^0 



-\/2Ar 



jliPjrr)-V2A, 



+4pl ( / dru(r) [(1 + x/2Ar) + V C^^^ (\/2Ar)"] ^ 
+ (Fo' + - Ml,)^j^ dru{r)j^{p^r)e-'''-'^) - ^pl[ £ dru{r){\ + MK*r)^^^^e'"''<'^) 



(A.6) 



with C^^^ = (Cf ^ - C^°^), C|^^ = (C^°^ - 2C;^'^^) and Cf ^ = \ Setting u(r) ^ r again that characterizes the weak 
binding limit one confirms by direct calculation that now 



Ra{p.) = ( 



2A2 - M]^, \ 4 



2A2+p2 



(A.7) 



instead of unity, where Pi(p^) is the form factor given by Eq. ^ evaluated at Mi = Mk* and = P^- Fo^' ^ — 420 
MeV the actual values are i?A = 0.39 versus Ff — 0.34 so the wave function at contact approximation continues to 
hold as well. Thus, depending on the value of the cutoff chosen, including a form factor in the transition amplitudes 
can reduce both the tttt and tt'^t] widths in Table [T] by an order of magnitude. 
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TABLE I: Calculated values (in MeV) of the isoscalar (eo) and isovector (ei) binding energies and resulting masses Mi — 

2Mk — ei and annihilation widths for /o(980) and ao(980) for varying cutoff A as given by the non-local meson exchange 
potential model. Theoretical and experimental values from various sources have also been listed. 

A eo Mq - -ir^^ ei Ml - |r^n„ 

415 21.3 971 - 54 i 15.5 976 - 13.7 i 

420 10.6 981 - 33 i 7.55 984 - 7.30 i 

425 4.03 988 - 18 i 2.56 989 - 3.74 i 

Local potential, Ref. f5] 18.6 981 — 25 i unbound - 

Linear cr model, Ref. [12| - 940 - 27 i - 983.4 - 69 i 

Fermilab E 791 Collaboration [23] - (975 ± 3) - (22 ± 2) i 

KLOE Collaboration [24] - - - (984.8 ± 1.2) - 61 i 

OBELIX Collaboration [25] - - - (998 ± 10) - 23 i 

Particle Data Group fll] - (980 ± 10) - (20 to 50) i - (984.7 ± 1.2) - (25 to 50) i 

TABLE II: Bound state properties of the isoscalar potential: depth Vo, Eq. (|15|l . effective mass ratio, Eq. (|13p . binding energy 
eo, wave function at contact amplitude 1/^(0), and the factor R of Eq. (|25p . 

A (MeV) V;)(0) (GeV) M^{0)/Mk eg (MeV) 7/>(0)/a4''^ R_ 

415 -0.23 2.18 21.3 0.126 1.041 

420 -0.21 1.93 10.6 0.096 1.038 

425 -0.19 1.75 4.03 0.072 1.034 
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FIG. 1: Critical values of 7(0), the inverse effective mass ratio versus the cutoff A for isospin 7 = and 7 = 1. Both curves 
remain positive in the common interval ^ 0.39 ^ A < 1.28 GeV. 
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FIG. 2: The direct and exchange diagrams that contribute to the additional four-point vertex, or box diagram f 6oa;(p', p; T'o), 
in the BS equation for KK scattering that arises from t-channel K* (892) exchange between the kaons, leading to tttt isoscalar 
or TT^r? isovector intermediate states. In the latter case only the first diagram is present. 
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FIG. 3: (a) Meson exchange diagrams for s-channel contributions to the Bethe-Salpeter transition ampHtude in Eq. Q. The 
thin and thick soHd lines indicate bare versus dressed intermediate mesons propagating between interaction vertices, shown 
as filled circles, (b) The meson proper selfenergy E. Only the tttt or tt^t] polarization loop contributions to E that include 
imaginary parts for dressing the bare a or ao meson propagators at s = Pg 



: 4Af|- are shown. 
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FIG. 4; The tttt inelasticity versus the total center of mass (cm.) energy. The filled circles with error bars show values of 
rjo{-KTT) taken from B. Hyams et al. [4ll ]. while the upper curve (a) gives the fit using the _fsr-matrix from [i^ . The lowest curve 
(c) shows the calculated inelasticity based on the effective range parameters of Eq. (|41|l . Curve (b) gives the result using the 
scattering parameters for the local potential in 
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